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ABSTRACT 

We present new evolutionary models for Type la supernova (SN la) progenitors, introducing mass-stripping 
effect on a main-sequence (MS) or slightly evolved companion star by winds from a mass-accreting white 
dwarf (WD). The mass-stripping attenuates the rate of mass transfer from the companion to the WD. As a 
result, quite a massive MS companion can avoid forming a common envelope and increase the WD mass up 
to the SN la explosion. Including the mass-stripping effect, we follow binary evolutions of various WD -i- MS 
systems and obtain the parameter region in the initial donor mass - orbital period plane where SNe la occur 
The newly obtained SN la region extends to donor masses of 6-7 Mq, although its extension depends on 
the efficiency of mass-stripping effect. The stripped matter would mainly be distributed on the orbital plane 
and form very massive circumstellar matter (CSM) around the SN la progenitor. It can explain massive CSM 
around SNe Ia/IIn(IIa) 2002ic and 2005gj as well as tenuous CSM around normal SN la 2006X. Our new 
model suggests the presence of very young (< 10** yr) populations of SNe la, being consistent with recent 
observational indications of young population SNe la. 

Subject headings: binaries: close — circumstellar matter — stars: winds, outflows — supernovae: individual 
(SN2002ic, SN 2005gj, SN 2006X) 



1. INTRODUCTION 

The nature of Type la su pernova (SN la) progenitors has 
not been clarified y et (e.g., lNiemever '& Hillebrandt 200^ 
iNomoto et al. II2OO0I) . although it has been commonly agreed 
that the exploding star is a mass-accreting carbon-oxygen 
white dwarf (Ch-O WD). For the exploding WD itself, the 
observed features of SNe la are better explained by the 
Chandrasek har mas s model than the sub-Chandrasekhar mass 
model (e.g.. lLivio I I2OOO). However, there has been no clear 
observational indication as to how the WD mass gets close 
enough to the Chandrasekhar mass for carbon ignition; i.e., 
whether the WD accretes H/He-rich matter from its binary 
companion [single degenerate (SD) scenario], or two C-i-O 
WDs merge [double degenerate (DD) scenario]. 

Recently, the following two important findings have been 
reported in relation to the SN la progenitors: (1) circumstellar 
matter (CSM) around the progenitors, and (2) a very young 
(< 10^ yr) population of the progenitors. 

Circumstellar Matter: In the SD scenario, H/He-rich 
CSM is expected to exist around SNe la as a result of 
mass transfer from the companion as well as the WD winds 
(e.g..l Nomotoll982HHachisu. Kato. & Nomotol 19 99a). Thus 
searching for H/He-rich CSM is one of the key o bservations 
to identify the progenitors ('e.g.. lLundqvist et al. ir2003i) . Re- 
cently detections of such CSM have been reported for several 
SNe la, i.e., observations of narrow H-emission lines in SNe 
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2002ic dHarnuv et al. ' '2003*) and 2005gi (Alderlin g etaL 
20061 iPrietoet al. . 2007) (Type la/II n or Ila (.Deng et al 



2004 ). fliermal X-rays fr om 2005ke dlmmler et al. I l2006h 



and Na I D lines in 2006X (iPatat et al. \\2001?it . 

The identification of SN 2002ic as an SN la has been 
confirmed by the recer it spectral comp arison between SN 
2005gj and SNe la (IPrietoetal. I 12007.) being against 
the Type Ic suggestion by iBenetti et aLl (l2006l) . Sev- 
eral CSM inte raction m odels suggested a 1-2 Mq CSM 
dChugai et al. 1 12004; No moto etaLll2005h . The evolution- 
ary origin of such a massive CSM has been explored by 
iLivio & Reiss I ((2003) based on a common envelope evolution 
model, by iHan & Podsiadlowski I (l2006l) from the delayed 
dynamical instability mod e l of b inary mass transfer, and by 
IWood-Vasey & Sokoloski I (l2006l) based on a recurrent nova 
model with a red giant companion. 

For normal SNe la, non-detection of radio has put the up- 
per limit of mass loss rate as M/vio < 1O~^M0 yr"', where 
v'lo = v/10 km s"' (Panagi a et al. 120 06^. However, the optical 
observations of SN 2006X have detected variable Na I D lines 
from CSM, whose expansion velocity an d mass have been es - 
timated to be vio - 10 and lO""* M© (iPatat et al. Il2007al) . 
Patat et al. have suggested that the CSM in SN 2006X origi- 
nated from the red-giant companion because of relatively low 
velocities. Comparing the SN 20 06X li ght curves with the 
other normal SNe la light curves, iWang et al.1 (l2007a ) sug- 
gested that the obvious deviation, the decline rate is slow- 
ing down in a later phase, can be explained by an interac- 
tion between ejecta and CSM or a light echo of circumstel- 



2 



Hachisu et al. 



lar/interstellar matter (see alsolWang et al. Il 2007bl). 
Young Population: According to Mannu cci et al. I (120061) . 

the present observational data of SNe la are best matched by a 
bimodal population of the progenitors, in which about 50 per- 
cent of SNe la explode soon after their stellar birth at the delay 
time of fdeiay ~ 10'^ yr, while the remaining 50 percent have a 
much wider distributi on of the delay time of fdeiay ^ 3 Gyr 
lAubourg et al. I (l2007l ) recently reported evidence for a short 
(less than 70 Myr) delay time component in the SN la popu- 
lation. In this paper, we define the term delay time as the age 
of a binary system at the SN la explosion , in order to compare 
our results with the earlier results (e.g., iGreggio & Renzini I 
[T98llGreggio]l200llMannucci et al. \\2(Mi . 

This kind of short delay times (fdeiay ^10^ yr) of SNe la 
have been suggested from the distribution o f SNe la relative 
to spi ral arms (e.g., B artun ov et al . 19 94; de lla Valle & Liviol 
11994 . Recently, Di Stef ano & Kong I (l2003l) reported, based 
on the Chandra data from four external galaxies: an elliptical 
galaxy (NGC 4967), two face-on spiral galaxies (MlOl and 
M83), and an interacting galaxy (M51), that in every galaxy 
there are at least several hundred luminous supersoft X-ray 
sources (SSXSs) with a luminosity of > lO-'^ erg s~' and that, 
in the spiral galaxies MlOl, M83, and M51, SSXSs appear 
to be associated with the spiral arms. The latter may indicate 
that SSXSs are young systems, possibly younger than 10** yr, 
and has some close relation to the young population of SNe 
la. 

The SD scenario has ever not predicted such young pop- 
ulations of fdeiay ~ 10^ yr. Corresponding to, at least, the 
zero-age main-sequence (ZAMS) stars at mass 5 -6_M^ 
(see, eTg., Li & van den Heuvel 1997; Hachisu et al. Ill999bt 
iLanger etal. I l2000t iHan & Podsiadlowski I ilOO^ In the 
present paper, we propose a scenario for such a young SN 
la population by introducing mass-stripping effect into binary 
evolutions. Mass-accreting WDs blow optically thick winds 
when the mass transfer rate to the WD exceeds the critical 
rate of Mer ^ 1 x IO'^'Mq yr"' dHachisu et al. 1 11996). The 
WD wind collides with the secondary's surface and strips off 
matter. When the mass-stripping effect is efficient enough, the 
mass transfer rate to the WD is attenuated and the binary can 
avoid the formation of a common envelope even for a rather 
massive secondary. 

The mass-strippin g effect on a MS compa nion has been 
first introduced by iHachisu & Kato I (I2003bll3) . who ana- 
lyzed two quasi-periodic transient supersoft X-ray sources, 
RX J0513.9-6951 and V Sge: RX J0513 shows a quasi- 
periodic oscillation between optical high (^ 100- 120 days) 
and low (~ 40 day s) states with an amplitude of 1 mag 
dAlcock etal. Ill996h . RX J0513 is X-ray b right only during 
the optical low states (iReinsch et al. Il2000h . IHachisu & Kato I 
(|2003b) proposed a model that the mass transfer is modulated 
by the WD wind because the WD wind collides with the com- 
panion and strips off its surface and attenuates the mass trans- 
fer rate. When the mass transfer rate decreases below the crit- 
ical rate M^r, the WD wind stops and supersoft X-ray turns 
on. This corresponds to an optical low sate. Then the mass- 
transfer rate recovers because of no attenuation by WD winds 
and the WD blows winds again. X-ray turns off and an opti- 
cal high state resumes and the binary starts the next cycle of 
quasi-periodic oscillation. Such a self-sustained model natu- 
rally explains major characteristics of quasi-periodic high and 
low states and this success encourages us to adopt the same 
idea in the evolution scenario of supersoft X-ray sources and 
SN la progenitors. 



In the present paper, we show that this mass-stripping effect 
derives (1) formations of circumstellar matter (CSM) around 
SNe la and (2) a very young population of SNe la. We sum- 
marize our basic treatments of mass-stripping effect and bi- 
nary evolutions in §2, and then show our numerical results 
and their relations to a very young population of SNe la in 
§3. In §4 we present the origin of CSM around SNe la based 
on our results and show a relation between the very young 
population of SNe la and their massive CSM. Discussion and 
concluding remarks follow in §§5 and 6. 

2. MASS-STRIPPING EFFECT AND BINARY EVOLUTION 

Strong winds from a mass-accreting WD collide with the 
companion star and strip off its surface. This mass-stripping 
effect plays an impo rtant role in binary evolutions (e.g. 
IHachisu etanil999al) . Here we reformulate its treatment in 
our binary evolution calculation. 

2.1. New Aspects of Binary Evolutions 

First we briefly introduce a new binary evolutionary process 
through stages (a)-(d) below (also shown in Fig. [TH-d), where 
(c) and (d) are new stages introduced by mass-stripping. 

(a) The more massive (primary) component of a binary 
evolves to a red giant star (with a helium core) or an AGB 
star (with a Ch-O core) and fills its Roche lobe. Mass transfer 
from the primary to the secondary begins and a common en- 
velope is formed. After the first common envelope evolution, 
the separation shrinks and the primary component becomes a 
helium star or a Ch-O WD. The helium star evolves to a Ch-O 
WD after a large part of helium is exhausted by core-helium- 
burning. We eventually have a close pair of a C+O WD and a 
main-sequence (MS) star as shown in Figure[T^. 

(b) After the secondary evolves to fill its Roche lobe, the 
mass transfer to the WD begins. This mass transfer occurs 
in a thermal timescale because the secondary mass is more 
massive than the WD. The mass transfer rate exceeds the crit- 
ical rate for the optically th ick wind to blow from the WD 
(IHachisu et al. l[T996l IT999a b). 

(c) Optically thick winds from the WD collide with 
the secondary surface and strips off its surface layer 
(IHachisu & Kato ll2003am . This mass-strippmg attenuates 
the rate of mass transfer from the secondary to the WD, thus 
preventing the formation of a common envelope for a more 
massive secondary in the case with than in the case without 
this effect. Thus the mass-stripping effect widens the donor 
mass range of SN la progenitors (see Fig. [3]below). 

(d) Such stripped-off matter forms a massive circumstellar 
torus on the orbital plane, which may be gradually expand- 
ing with an outward velocity of ^ 10- 100 km s"' (Fig. [TJl), 
because the escape velocity from the secondary surface to L3 
point is v'esc ^ [(0L3-0Ms)GM/fl]'/^ ^ 100 km s"' (see be- 
low). Subsequent interaction between the fast wind from the 
WD and the very slowly expanding circumbinary torus forms 
an hourglass structure (Fig. [T];-d). 

2.2. Formulation of Mass-stripping 

Fast strong winds collide with the companion as illustrated 
in Figure[T] The companion's surface gas is shock-heated and 
ablated in the wind. We estimate the shock-heating by assum- 
ing that the velocity component normal to the companion's 
surface is dissipated by the shock and the kinetic energy is 
converted into the thermal energy of the surface layer. The 
heated surface layer expands to be ablated in the wind. 
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Fig. 1. — A schematic configuration of a binary evolution including mass-stripping effect, (a) Here we start a pair of a C+O WD and a more massive 
main-sequence (MS) star with a separation of several to a few tens of solar radii, (b) When the secondary evolves to fill its Roche lobe, mass transfer onto 
the WD begins. The mass transfer rate exceeds a critical rate for optically thick winds. Strong winds blow from the WD. (c) The hot wind from the WD hits 
the secondary and strips off its surface, (d) Such stiipped-off material forms a massive circumstellar disk or torus and it gradually expands with an outward 
velocity of 10— 100 km s"'. The interaction between the WD wind and the circumstellar torus forms an hourglass structure. The WD mass increases up to 
Mia = 1-38 Mq and explodes as an SN la. When we observe the SN la from a high inclination angle such as denoted by "line of sight," circumstellar matter 
(CSM) can be detected as absorption lines like in SN 2006X. 

To obtain the mass stripping rate, we use th e same formula- 
tion proposed by iHachisu & Katol (12003 b'' We equate the 
stripping rate times the gravitational potential at the compan- 
ion surface to the net rate of energy dissipation by the shock 
as: 



effect of Roche lobe overflow from the L3 point. Then the 
stripping rate is estimated as 



^^.strip =ClMwind, 



GM -12 

((PL3-(PMS)-Msmp= -V ■ Vefl ■ g(l) ■ M„ind, (1) 

a 2 
where M = Mwd+A^ms, A^wd is the WD mass, Mms is the 
main-sequence companion mass, a is the separation of the bi- 
nary; (pMS and 0L3 denote the Roche potential (normalized by 
GM /a) at the MS surface and the L3 point near the MS com- 
panion, respectively; v is the WD wind velocity, yyeff is the ef- 
ficiency of conversion from kinetic energy to thermal energy 
by the shock, g{q) is the geometrical factor of the MS surface 
hit by the wind including the inclination (oblique shock) ef- 
fect of the wind velocity against the companion's surface (see 
IHachisu. Kato. & Nomotol il999ai for more details on g{q)), 
and q = M2/M1 = Mms /-^wd is the m a ss ratio . Here we mod- 
ified equation (21) of IHachisu et al. I (Il999al) to include the 



where 



ri^n ■ gig) f v^a \ 



(2) 
(3) 



Here we assume that the WD wind is spherically symmetric. 
If the asphericity of the WD wind is not so large, having a lat- 
itudinal (6'-angle) dependency like a broad angle jet, we have 
a different form of g{q) and its value may be much smaller 
than that for the spherically symmetric WD winds. We also 
assume T^etf = 1 in the present calculation. When the wind ve- 
lo city is as fast as 4,0 00 km s~', we have c\ ~ 10 as estimated 
by IHachisu & Katol (2003b). Although ther e is a large ambi- 
guity in this kind of parameterization as ci . IHachisu & Katol 
t2003biCi) found the best fit models with c\ = 1.5-10 for 
RX J0513.9-6953 and cj = 7-8 for V Sge. We thus assume 
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Fig. 2. — SN la evolutions for two typical cases of WIND and CALM, (a) Case WIND: starting from Mwd.o = 10 Mq, M2.0 = 5.0 Mq, and Pq = 2.15 days 
with ci = 3, the WD reaches the SN la explosion in the wind phase at / = 6.57 X 10^ yr. The WD mass (Mwd), secondary mass (M2), mass loss rate from the 
secondary (M2), WD wind mass loss rate (M^inj), radius of the secondary (R2), effective radius of the Roche lobe for the secondary (SJ), and orbital period 
(Porb) are plotted. Only the orbital period is multiplied by four to easily see its change, (b) Case CALM: starting from Mwd.o = 1-0 Mq, M2,o = 5.0 Mq, and 
Po = 6.79 days with ci = 3, the WD reaches the SN la explosion but in an SSXS phase without winds at f = 6.93 X 10^ yr. The WD wind stops at f = 5.5 X lO' yr 
Even after that, the WD loses its mass due to weak helium shell flashes IKato & Hachisu 1999). Here M^ind includes an average mass loss rate by helium shell 
flashes and thus does not become zero after the optically thick wind of steady hydrogen shell burning stops. Values of the secondary radius (R2) and the Roche 
lobe radius for the secondary (RJ ) are divided by two to squeeze them into the figure. 



ci = 1, 3, and 10 to examine the dependence on the mass- 
stripping effect, because the essential ambiguity of our for- 
mulation is included in the ci parameter. 

When winds blow from the WD and strip off the compan- 
ion's surface, the change of the separation, a, is calculated 
from 



d Ml +M2 Ml Ml J 

- = — --2 — --2 — - + 2- 

a M1+M2 Ml M2 J 

Ml +M2 ^Mi ^Mi 

M1+M2 Ml M2 

M1+M2 • „ • s 



mated as 



/ = 



1 



1+q 



(5) 



(6) 



(4) 



where Mi = Mwd, M2 = Mms, and are the specific an- 
gular momenta of the WD wind and the stripped-off matter, 
respectively, in units of a^ilorb with fiorb being the orbital an- 
gular velocity. Since the WD wind is much faster than the 
orbital motion, the wind cannot get angular momentum from 
the orbital torque during its journey, so that the wind has the 
same specific angular momentum as the WD, which is esti- 



The ablated gas from the companion is assumed to have the 
angular momentum at the companion's surface. Then we have 
a numerical factor of 

which was given in Table 1 of iHachisu et al. I (| 1999a) and is 
rather small compared with (See IHachisu et al71ll999ai 
for more details of l^.) 

2.3. Modified Mass Transfer Rate 

We have followed binary evolutions from the initial state 
of (Mi,o, M2.Q, Pq), i.e., (MwD.o, Mm&.q, Pq), where Po is the 
initial orbital period. Here, the subscript naught (0) denotes 
stage (a) in Figure [T] that is, before the mass transfer from 
the secondary starts. The radius, R2{M2,t), and luminosity, 
L2(M2,t), of stars which have slightly evolved off from the 
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zero-age main-sequ ence (ZAMS), are c alculated using the an- 
alytic form given by lTout et al.1(ll997l) . 

The mass transfer proceeds on a thermal time scale when 
the mass ratio M2/M1 exceeds 0.79. We approximate the mass 
transfer rate as 



-Mo 



M2 
tkh 



■ max 



Crl ~ Cms 
Cms 



,1 



where trh is the Kelvin-Helmholtz timescale given by 



Tkh « 3 X 10^ yr 



(7) 



(8) 



(e.g..' Paczvnskilll971h . and Crl = d\ogR*/d\ogM and Cms = 
dlogRMs/dlogM are the mass-radius exponents of the inner 
critical Roc he lobe and the main seq uence component, respec- 
tively (e.g., |HiinmingA^bbink][i987i). The effective ra- 
dius of the inner critical Roche lobe, R*, is calculated from 
Eggleton's (1983) empirical formula, i.e., 

a 0.6^2/3 + in(i+^i/3)' 

where ^ = M2 /Ml. 

When the mass transfer rate to the WD exceeds a critical 
value, which is given by 

'MwD 



iO.75 X 10" 



Mr 



0.4 Mq yr" 



(10) 



for the solar composition (hydrogen content of X = 0.7 and 
metallicity of Z = 0.02), the WD blows a wind with a mass 
loss rate of Mwind (< 0). This critical rate of M^r is the same 
as the critical rate for mass-accreting WDs to expand to a giant 
size, i.e., Mrg (see Nomoto et al.1l2007l for the recent calcu- 
lation of Mrg). The mass loss from the WD also occurs during 
the hydrogen shell flashes when -M2 < Mstabie, where Mstabie 
is the lowest rate for steady hydrogen burning and given by 
equation 



Mutable « 0.31 X 10" 



M 



WD 



Mfr 



0.54)MQyr-' (11) 



(|Nomoto et al. "2007). When Mstabie < -M2 < Ma, we have 
no mass loss associated with steady hydrogen shell-burning 
but have mass loss by helium shell flas hes. This mass 
loss play some role in the binary evolution ( Kato & Hachisu I 
119991) . Therefore, M^ind is the summation of the optically 
thick wind mass loss, hydrogen shell flashes, and hehum sheU 
flashes. 
We have the relation 

Mi +M2= Mwind +Mstrip, (12) 

from the total mass conservation, thus defining the net mass 
transfer rate to the WD as 



M„ 



transfer 



Mstrip-M2=Mi-M, 



wind ; 



(13) 



where signs of Mnansfei > 0, Mstrip < 0, M2 < 0, Mi > 0, and 
Mwind < should be noted. If M2 is given, we have the net 
mass transfer rate of 

,v /(ciMci-M2)/(ci + l), for -M2 >Mci- 

^"=1 -M2, for-M2<M„' ^^^^ 

where we use equations dU, ( fTST l, and a relation of 

-Mwind = Mtransfer-Mcr, (15) 



for -M2 > McT- Other treatm ents for binary evolutio n are es- 
sentially the same as those in lHachisu et al. I (Il999bh . 

Figure |2] shows two typical evolutionary sequences that 
demonstrate the effects by the modified mass transfer rate, 
M2, in equation 

(a) Starting from Mwd.o =1-0 Mq, M2.0 = 5.0 M©, and 
P() = 2.15 days with ci = 3, the WD reaches the SN la ex- 
plosion in the wind phase (Case WIND) at f = 6.57 x 10^ yr 
after the secondary fills its Roche lobe. The WD increases 
its mass (Mwd) up to Mia = 1-38 Mq to explode as an SN la. 
The secondary mass (M2) decreases to 2.01 Mq at the explo- 
sion. Both the mass decreasing rate of the secondary (dashed 
line labeled M2) and the WD wind mass loss rate (dashed line 
labeled Mwind) are also decreasing rapidly especially in the 
early phase of f < 1 x 10^ yr. This is because -M2 is large 
and the mass transfer rate, Mtransfer, is large during this phase, 
and as a result, both the WD wind mass loss rate, Mwind, and 
the stripping rate, Msmp, are also large. Shortly after this early 
phase, the Roche lobe's mass-radius exponent, Crl, becomes 
smaller than the secondary's mass-radius exponent. Cms, that 
is, Crl "Cms < 0. This gives -M2 = M2/TKH from equation 
(Q. We keep this mass transfer rate as long as the secondary 
overfills the Roche lobe, i.e., R2 > /?2- Figure|2^, we plot 
the secondary radius (the red line labeled R2) and the Roche 
lobe radius for the secondary component (the blue line labeled 
R'2) to show the condition of R2 > R2 during the evolution. 

(b) Starting from Mwd.o = 1-0 Mq, M2,o = 5.0 Mq, and 
Po = 6.79 days with ci = 3, the WD reaches the SN la explo- 
sion but in a phase of no winds (Case CALM) at f = 6.93 x 10^ 
yr after the secondary fills its Roche lobe. In this case the evo- 
lution of the mass transfer rate is different from Case WIND 
above. With -M2 =M2/tkh for Crl < Cms in equation (|7]), the 
secondary eventually underfills the Roche lobe, i.e., R2 < /?2- 
This can be seen in Figure |2j5, where the line of R2 crosses 
the line of at f ^ 1 x 10^ yr. This is because the stripped 
matter has rather low specific angular momentum (eq. ||6l), 
so that the binary separation hardly shrinks or even increases 
as seen from the temporal increase in the orbital period in 
Figure In realistic binary evolutions, the mass transfer 
is tuned in a way that the secondary radius is always equal 
to the Roche lobe radius for the secondary, i.e., R2 = 
Therefore -M2 is drastically decreased after f ^ 1 x 10^ yr, 
as shown in Figure |2j3. Thus, the optically thick WD wind 
stops at f = 5.5 x 10^ yr. In such a low mass transfer phase as 
Mtransfer ~ 1 X lO'^M© yr~'. Weak helium shell flashes occur 
and play an important role as a mass loss mechanism. This 
helium flash wind also strips off the secondary surface, thus 
working as a stripping effect. We introduce mass-stripping 
effect by these helium shell flashes into our binary evolution. 
Very small but finite Mwind in Figure|2]3 (after winds stop) rep- 
resents the mass loss from the WD at helium shell flashes and 
M2 includes the ensuing mass-stripping from the secondary. 

3. YOUNG POPULATION TYPE lA SUPERNOVAE 

Based on the bin ary evolution scenario proposed by 
iHachisu et al. I (Il999allbl) . we have followed binary evolutions 
starting from stage (b) in Figure[Tl that is, just when the com- 
panion evolves to fill its Roche lobe. The main difference 
from the previous work cited above is the inclusion of mass- 
stripping effect. Our results are shown in Figures [3l-[T0l 

Figure [3] shows the parameter regions that produce SNe 
la (SN la region) in the logfo-M2,o (the initial orbital pe- 
riod and the initial secondary mass) plane for the WD + MS 
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Fig. 3. — The initial parameter regions producing SNe la are plotted in the logPo~^2,o (orbital period — donor mass) plane for the WD + MS systems 
with various mass-stripping factors, cj . Thick solid: cj = 10. Medium solid: ci = 3. Thin solid: C| = 1. Dotted: cj = 0. The (red) hatched region indicates a 
region with a short delay time (/delay < 100 Myr) for the case of cj = 10. The region extends to the more massive donors for the larger ci . Tw o supersoft X-ray 
sources, RX J05 13.9—6951 (open circle) and V Sge (filled circle), are plotted, mas ses of which are esti mated to be 2.7 Mrr, j Hachisu & Kato I 2003h) and 3.5 Mp , 
(Hachisu & Kato 2003c), orbital periods of which are determined to be 0.76 days jPakull et al. 119931) and 0.51 days Merbig et al. ll965tlPatterson et al. 11993) . 
respectively. The position of V Sge suggests that C[ > 0. 
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TABLE 1 

Three typical cases of SN Ia explosion 



case 


wind 


H burning 


CSM 


pre-SN history 


SN la 


delay time 


immediate radio/X-ray 


WIND 


wind 


steady 


massive: near 


WIND (V Sge type) 


Ila (02ic-like) 


young 


yes 


CALM 


no wind 


steady 


thin: far 


WIND^SSXS 


normal la 


young 


no (~ 10-lOOyr) 


RN 


no wind 


flash 


very thin: many shells 


WIND^SSXS^RN 


normal la 


broad 


no(~ 100- 1000 yr) 










or SSXS^RN 
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Fig . 5 . — The pai'ameter region that produces SNe la is plotted in the log P- (orbital period — donor mass) plane for the WD + MS system. Here we assume 
^WD.o = 1-1 Mq for the initial white dwarf mass. The initial WD + MS system inside the region encircled by the (red) thin soUd line (labeled "initial") will 
increase its white dwarf mass up to the critical mass (Mi^ = I.SSMq) for the SN la explosion to occur. The final state of the WD + MS system in the logP-Mj 
plane just before the SN la explosion is encircled by the (blue) thick solid line (labeled "final"). The final state of the WD just before the SN la explosion is 
specified by one of wind (open circle), steady H-burning (filled triangle), or recurrent nova (open square) phase. An hatched region indicates a region in which 
the progenitor explodes in a delay time of ?delay < 100 Myr. Dashed line: in a delay time of 200 Myr. Dotted line: in a delay time of 400 Myr. Currently known 
positions of three recuirent novae are indicated by a star mark (★) for U S co (e.g., Schaefer & Ringwald, 1995; Hachisu et al. 2000a b), and by aiTows for the 
other two recuirent novae, V394 CrA tSchaefer. 19901) and CI Aql (Mennickent & Honevcutt .1995i) , of unknown companion masses. The WD ma sses of U Sco 
and V394 CrA were estimated to be 1.37 Mm Machisu et al. r2 000a; Hachisu & Kato 2000) while that of CI Aql was 1.2 jHachisu & Kato 120 03a). 
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Fig. 8. — Same as Fig. |6] but for an initial white dwarf mass of Mwd.o = 1.0 Mq. Large difference in the border of WIND-CALM comes from the fact that 
the secondary considerably overfills the Roche lobe, i.e., i?2 > at the SN la explosion in the Case WIND. 



system. Here the initial white dwarf mass is assumed to be 
My^Dfi =1.0 Mq. The white dwarfs within these SN la re- 
gions will increase their mass, Mwd, up to the critical mass 
(Mia =1.38 Mq) for the SN la explosion to occur. 

The SN la region in the logPo-A^2.o plane is enclosed by 
four boundaries. (1) The left boundary is given by the mass- 
radius relation for the zero-age main-sequence stars. (2) The 
lower boundary is set by strong nova explosions, below which 
A^transfer ^ 1 X 10~^Mq yr~' and the resultant nova explosion 
ejects most of the accreted matter, thus preventing the WD 



mass from increasing. (3) The upper boundary is limited by 
the formation of a common envelope. Here we assume that 
a common envelope is formed when Mtransfer ^ 1 x 1Q~'^Mq 
yr~' because R \.r\h ^,a^ 10 Rq for such a high Mtransfer (see 
iHachisu et al. |[l9 99b, for more details). (4) The right bound- 
ary corresponds to the end of central hydrogen burning of the 
MS companion: after that, it shrinks and underfills its Roche 
lobe. 

In Figure[3] the SN la regions for the various mass-stripping 
factor ci = 10, 3, and 1 are encircled by the thick, medium, and 
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Fig. 10. — Same as Fig.|6] but for an initial white dwarf mass of Mwd.o = 0.9 Mq. 



thin solid lines, respectively, and the no stripping case (ci = 0) 
by the dotted line. The position of the Galactic supersoft X- 
ray source V Sge is clearly outside the SN la region for ci = 0, 
but inside the SN la region if ci > 0. For larger ci, the SN 
la region extends to more massive M2.0, because the stronger 
mass-stripping leads to the lower mass transfer rate, Mtmnsfer, 
from the secondary to the WD (see eq. lfl4l '). thus prevent- 
ing the formation of a common envelope for larger M2,o- As 
shown in this figure quite massive secondaries produce SNe la 
(e.g., M2fl = 7.5 Mq for ci = 10) for the strong mass-stripping 



case of ci > 3. 

Such WD -I- MS systems with a massive MS secondary con- 
sist of a very young population of SNe la. We show the region 
of short delay times, fdeiay < 100 Myr, by the red shadow in 
Figures[3]|ll|5]|2l and|9] Figure |4] shows the SN la regions for 
different initial WD masses, Mwd.o = 0.7, 0.8, 0.9, 1.0, and 
1 . 1 Mq . The red (sparse) and blue (dense) hatched regions in- 
dicate the delay time of fdeiay < 100 Myr for Mwd,o = 1-1 Mq 
and 0.7 Mq, respectively. 

We apply the present result to equation (1) of 
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Ilben & Tutukovl ( 11984 . i.e., 

iy = Q.2-Aq-J^ _.Alogfl yr"', (16) 

where A^, A log a, Mi, and M„ are the appropriate ranges 
of the mass ratio and the initial separation, and the lower and 
upper limits of the primary mass for SN la explosions in solar 
mass units, respectively. We then estimate the SN la birth rate 
in our Galaxy as i^wD+MS ^ 0.004 yr~ ', which is consistent 
with the observation (Cappell aro et al. I[l9 99). 

On the other hand. Hachisu et al. I(ll999al) proposed another 
channel to SNe la, the symbiotic channel, binary of which 
consists of a white dwarf and a red giant (WD + RG), and 
estimated its birth rate to be t'wo+RG ^ 0.002 yr"'. 

Assuming the initial distribution of binaries given by 
equation (fT6] l at the burst of star formation (single event), 
we estimate the delay time distribution of SNe la for the 
WD + MS systems in Figure [TT] The number ratio of 
these young populations is calculated for 10 bins of delay 
time, (0.025,0.05), (0.05,0.1), (0.1,0.2), (0.2,0.4), (0.4,0.8), 
(0.8,1.6), (1.6,3.2), (3.2,6.4), (6.4,12.8), and (12.8,25.6) 
Gyr. The number ratio with f delay < 100 Myr and f delay < 
200 Myr are about 50% and 80%, respectively, of the to- 
tal SNe la coming from the WD + MS system, which is 
consistent with the recent observational s uggestions (e.g., 
iMannucci et al. ll2006l:lAubourg et aini2007l) . 

Short delay times (fdeiay ^10^ yr) of some SNe la have been 
suggested from the distribution of SNe la relative to spira l 
arms (e.g.. .Bartunov et al. Ill994t Idella Valle & Livio 111994 . 
iPetrosian et al. I (l2005h reported that about 30-40% of SNe la 
are associated with spiral arms in their samples, being consist- 
ing with our results. Mannucc i et al. I (j2006) have suggested 
that the delay time distribution function of SNe la has a bi- 
modality, one for young population (fdeiay 100 Myr) and the 
other with a broad distribution over 3 Gyr Our delay time 
distribution function has a peak around fdeiay < 100 Myr from 
the WD -I- MS systems and a broad d istribution from the WD 
-I- RG systems (iHachisu et al. lfT999al) as shown in Figure [12] 

4. FINAL STAGE OF BINARY EVOLUTION AND CIRCUMSTELLAR 

MATTER 

The final state of the WD depends mainly on the mass trans- 
fer rate Mtransfer f rom the dono r star to th e WD at the S N la ex- 
plosio n dNomoto 1982; Hachisu et al. iri999a; Nomoto et al. I 
l2007h . As shown in Figure |2] M2 drops quickly in the early 
stage and then slows down to almost a constant value. At 
least, in the early phase, the mass transfer proceeds on a ther- 
mal time scale, represented by the second term of equation 
(|7]l, when the mass ratio M2/M1 exceeds 0.79. So we approx- 
imate the mass transfer rate as 

-M2 ~ ~ 3 X 10 Mq yr 

(17) 

By applying the approximate M2-L2 relation of L2 oc M™ , 
where m ^ 4 for the 1.5-3 Mq zero-age main-sequence 
(ZAMS) stars or m - 3 .5 for the 3 - 7 Mq ZAMS stars, 

-M2 0c/?2Mr^ (18) 

Thus -M2 decrease as M2 decreases. 

Figures ISHTOl show the SN la regions in the logf -Md (or- 
bital period — donor mass) plane for the initial WD -1- MS 
system (encircled by the red thin line and labeled "initial") 



as well as the final state at the SN la explosion (encircled 
by the blue thick line and labeled "final"). Here we assume 
ci = 3 and Mwd.o = 11, 1-0, and 0.9 Mq. In these figures, we 
distinguish three final states just before the SN la explosion, 
i.e., optically thick WD wind phase (WIND: open circles), 
steady hydrogen burning phase without optically thick winds 
from WDs (CALM: filled triangles), and recurrent nova (RN) 
phase (RN: open squares). The characteristic properties for 
these three progenitor stages are summarized in Table [T] and 
the corresponding binary parameters are tabulated in Table |2] 

4.1. Case WIND 

When the mass transfer rate from the secondary continu- 
ously exceeds the critical rate of equation (fTol i until the fi- 
nal stage, the WDs explode during the wind phase (Fig. |2^). 
Therefore, we call this Case WIND. Case WIND is real- 
ized in the region of M2,o > 3 Mq and P2.0 < 2 days for 
MwD.o =1.1 Mq and 1 .0 Mq (open circle), but no Case WIND 
exists forMwD.o < 0.9 as shown in Figures|5]|7] and|9] 

The stripped-off matter from the companion can easily 
amount to AMstrfp ^ 1-2 Mq and even reach 3-4 Mq as seen 
from the donor mass difference AM2 between the "initial" and 
the "final" in Figures |5]|7] and|9] More precisely, AM2 con- 
sists of three parts, the stripped-off mass AMstdp, the accreted 
mass by the WD AMi, and the mass ejected by the WD wind 
AMwind, i.e., M2= Msmp+Mwind-Mi from equation ( fT2b . This 
can be approximated as M2 ~ Mstrip +A^wind = ( 1 + 1 /ci )Mstrip = 
4/3Mstrip because Mi < -M2, so that AMsmp « 3/4 AM2 for 
ci = 3. 

The stripped-off material forms CSM very near the SN la. 
We expect that stripped-off matter did not go away from the 
system because the velocity of stripped-off matter may not 
exceed the escape velocity of the binary system. Then the 
SN la undergoes circumstellar interaction as observed in Type 
la /Iln (or 11a) SNe 2002ic and 2005gi. 

lAlderhng et al. I (l2006i) suggested that the host galaxy of 
SN 2005gj had a burst of star formation 200 ± 70 Myr ago. If 
the progenitor of SN 2005gj was born at that time, its delay 
time is consistent with our Case WIND as shown in Figure 

4.2. Case CALM 

When the mass transfer rate from the secondary is below 
the critical rate for optically thick winds but above the lowest 
rate of steady hydrogen burning, i.e., MstaWe < A/transfer < M^, 
the WDs undergo steady H-burning at the time of SN la ex- 
plosion (filled triangles in Figs. l5HT0b. We call this Case 
CALM because no optically thick winds occur. The WDs are 
observed as supersoft X-ray sources (SSXSs) until the SN la 
explosion. The stripped-off material forms CSM but it has 
been dispersed too far to be detected immediately after the 
SN la explosion. 

The CALM case is realized in the region of M2.0 > 3 Mq 
and P2.0 ^ 2 days for Mwd.o = 11 Mq and 1 .0 Mq in Figures 
|5] 111 and |9] where Mtransfer in the early phase is much larger 
than that of ^2,0 ^ 2 days, because in equations (|7]i and ( fTTb . 
R2 and L2 are much larger than those for ^2,0 ^ 2 days. Then 
Mtransfer is much larger, thus much more mass had been lost 
in the earlier phase. As a result, the wind phase finishes at 
an earlier time even for the same initial mass M2,o as seen in 
Figures|2^ and|2]3. Therefore, at the SN la explosion, no wind 
occurs. 
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TABLE 2 

Initial parameters for Three SN Ia explosions 



WD mass 


secondary mass 


orbital period 


case 


pre-SN history 


SN la 


(Mq) 


(Mq) 


(days) 












1.0-1.1 


3-6 


~ 0.5-2 


WIND 


WIND 






Ila (02ic-like) 


1.0-1.1 


3-6 


~2-10 


CALM 


WIND- 


►SSXS 




normal la 


1.0-1.1 


2.2-3 


~ 0.5-4 


CALM 


WIND- 


►SSXS 




normal la 


1.0-1.1 


1.8-2.2 


--0.5-2 


RN 


WIND- 


»SSXS- 


»RN 


normal la 










or SSXS 


-^RN 






0.9 


2.5-5 


~ 0.5-6 


CALM 


WIND- 


»SSXS 




normal la 


0.9 


2.0-2.5 


~ 0.5-2 


RN 


WIND- 


►SSXS- 


»RN 


normal la 


0.8 


4-5 


~ 1-3 


RN 


WIND- 


►SSXS- 


»RN 


normal la 


0.8 


4-5 


~0.5-l 


CALM 


WIND- 


►SSXS 




normal la 


0.8 


2.5-4 


~ 0.5-2 


RN 


WIND- 


►ssxs- 


»RN 


normal la 


0.7 


3-4.5 


~0.5-l 


RN 


WIND- 


►ssxs- 


»RN 


normal la 



In the region of M2X) < 3 Mq for Mwd,o =1-1 Mq and 
1 .0 Mq in Figures |5]-[8] (filled triangles or open squares), M2 
decreases to as small as the primary Mwd, i-e., the mass ra- 
tio of ^ ~ 1, at the SN la explosion, which corresponds to a 
lower part of the "final" region. Then the mass transfer rate 

decreases down to Mu-ansfer < M^r or even Mu-ansfer < A^stabie be- 
cause L2 is smaller for the smaller M2 even if the mass transfer 
itself is proceeding on a thermal time scale. The wind phase 
has ended earlier than for Mofi > 3 Mq. 

The border between Case WIND and Case CALM can 
be simply estimated from the condition Mnansfer = Mq and 
R2 = R2 at the SN la explosion. We have calculat ed M2 from 
eguation ([TtT i with R2 and L2 being taken from iTout et aO 
(|T997) (a single star evolution). This line agrees reasonably 
with the border of WIND-CALM in Figure but largely de- 
viates from it in Figure |8] This is because the secondary con- 
siderably overfills the Roche lobe, i.e., R2 > R^, at the SN la 
explosion for Mwd,o = 1-0 Mq. 

For MwD.o ^ 0.9 Mq (filled triangles or open squares), M2 
decreases at the SN la explosion as shown in Figure|9]("final" 
region) and Mtransfer decreases to be lower than M^ mainly 
because the time for the WD to reach Mia = 1 -38 Mq is longer 
and much more mass is lost during the evolution. The wind 
phase has ended before the SN la explosion. 

For a typical case of ci = 3, Mwd.o = 0.9 Mq, M2.0 = 4.0 Mq, 
and P() =1.3 days in Figure |9] the WD explodes as an SN 
la at f = 9 X 10^ yr after the secondary fills its Roche lobe. 
The wind has already stopped 3 x 10^ yr ago (duration of the 
WIND phase, Afwind = 6x10'' yr, and duration of the CALM 
phase, Afcaim = 3x10^ yr), so that the inner edge of stripped- 
off material has already gone to (10- 100 km s"') x (3 x 
10^ yr) ~ lO'"*- 10^0 cm from the SN la. Therefore, it takes 
about 10- 100 yr for the SN la ejecta to reach the inner edge 
of stripped-off matter. We do not expect radio or X-ray until, 
at least, 10-100 yr after the explosion. Thus the resultant 
SNe la are mostly "normal." The duration of CALM phase 
is typically a third or fourth of the total evolution time to the 
SN la. These long durations of optically thick wind phases 
may reduce the statistical number of luminous supersoft X- 
ray sources because the photospheric temperature of the WD 
is lower than ^ 10 eV and not luminous in supersoft X-ray. 

The decline of SN 2006X light curves is slowing down in 
a later phase compared with the other normal SNe la light 
curves, suggestin g an interaction b etween the ejecta and CSM 
in a later phase ( Wang et al. 2007a) or a light echo of circum- 
stellar/interstellar matter (Wang et al. 2007b) . This happens 
if SN 2006X is placed at the border between our Case WIND 
and Case CALM since the innermost part of slowly expand- 



ing circumstellar matter has not yet moved far away. A X-ray 
detected SN 2005ke may also belong to the same category 
dlmmler et al. II2006I) . 

For the progenitor of SN 2006X, iPatat et al. I (l2007al) sug- 
gested a WD + RG system like RS Oph from the circumstellar 
matter (CSM) absorption lines. Here we suggest that a WD + 
MS system (like U Sco) may better explain a continuous ve- 
locity distribution (from -30 to -150 km s"') of the CSM 
absorption lines by the stripped matter with continuous veloc- 
ity distribution (see Fig. [TJl). In this connection, very recent 
report of the Na I D circums tellar lines of RS Oph during the 
2006 outburst is suggestive (Ilijimall2007h . These lines indi- 
cate no continuous distribution as observed in SN 2006X but 
a narrow velocity component of -36 km s~' against RS Oph 
that is attributed to the red giant cool wind. 

Recently negative detections of time-variable Na I D lines 
have b een reporte d for two SNe la 2000cx jP atat et alj 
I2007b|) and 2007af (ISimon et al. Il2007h . iPatat et al. ■ (,2007bh 
and lSimon et al. I (l2007h suggested a possibility that the distri- 
bution of CSM is torus/disk-like as illustrated in Figure [T| In 
such a case, variable Na I D lines would not be observed if the 
line of sight is perpendicular to or off the orbital plane. Since 
the hot WD winds have a large velocity of > 1000 km s~\ the 
CSM formed by hot winds quickly diffuse away and is too 
tenuous to be detected; 

Recently, Badenes et al. I (l2007h reported that the fast WD 
wind of V > 200 km s"', which excavates its circumstellar 
medium and forms a large cavity around an SN la, is incom- 
patible with the X-ray emission from the shocked ejecta in our 
Galaxy (Kepler, Tycho, SN 1006), Large Magellanic Cloud 
(0509-67.5, 0519-69.0, N103B), and M31 (SN 1885). We can 
avoid this difficulty if the stripped-off matter has a velocity of 
10- 100 km s-i. 

4.3. CaseRN 

When the mass transfer rate from the secondary is below the 
lowest rate of steady hydrogen burning, i.e., Mtransfer < Mstabie, 
hydrogen shell burning is unstable to flash and recur many 
times in a short period as a recurrent nova (RN) (the open 
squares in Figs. ISHTOll. We call this Case RN. The recurrent 
nova U Sco, one of the candidates of SN la progenitors, is in 
the middle of the "final" region (Hachisu et al. 2000a b). The 
resultant explosions are "normal" SNe la. 

A simple estimation gives the border between Case CALM 
and Case RN, Mtransfer = Mstabie and R2 = R*2 at the SN la explo- 
sion. Here we cal culate M2 from equ ation (fTTT i with R2 and L2 
being taken from lTout et al. I d 19971) (a single star evolution). 
These lines agree reasonably with the border of CALM-RN 
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in FiguresillS] and\M 

For a typical Case RN of Mwd,o = 1 .0 Mq, M2,o = 2.0 Mq, 
f = 118 days with ci = 3, the WD undergoes the SN la ex- 
plosion in the recurrent nova phase at f = 9.49 x 10^ yr af- 
ter the secondary first fills its Roche lobe. The WD wind 
stops at f = 4 X 10^ yr and the stable hydrogen burning ends 
at f = 8.6 X 10^ yr. During the last 10^ yr in the recur- 
rent nova phase, the secondary loses ^ 0.022 Mq, of which 
the WD accretes 0.017 Mq. Therefore the stripped-off mat- 
ter in the recurrent nova phase is very small. On the other 
hand, the stripped-off matter in the early wind phase amounts 
AMstrip ~ 0.15 Mq, which has already been far from the SN 
at the SN la explosion, i.e., (10-100 km s"') x (5 x lO^* yr) 
= (1 - 10) X 10''' cm. It takes about 100-1000 yr for the SN 
ejecta to reach the stripped-off matter. These features are 
summarized in Table [T| 

5. DISCUSSION 

5.1. Mass-stripping Effect and Modulated Mass Transfer 

Rate 

As mentioned in §1, the existence of mass-stripping ef- 
fect has been demonstrated by Hachisu & Kato (2003b c). 
They analyzed two quasi-periodic transient supersoft X-ray 
sources, RX JOS 13. 9 -6951 in the Large Magellanic Cloud 
(LMC) and V Sge in our Galaxy. Especially V Sge shows 
the following key observational features: (1) V Sge exhibits 
long-term transitions between optical high (its brightness of 
y ^ 11 and its duration of ^ 180 days) and low (V ~ 12 
and 120 days) states with total durations of 300 days 
(see, e.g.. ISimon & Mattel 1 1 19991 for the long-term behav- 
ior). (2) Very soft but very weak X-rays are detected only in 
the lo ng-term optical low state (e.g., Greiner & van Teeseling 
Il998l) . (3) Radio observations indicate a wind m ass-loss rate 
as large as ^ 10"^M q yr"' (,Lockley. Eyres. & Woodlll997l: 



lLocklevetal~lll999l). 

iHachisu & Kato 1 (120033 ) explained these features based on 
the mass-stripping effect; the mass transfer to the WD is 
modulated by the WD wind because mass-stripping attenu- 
ates the mass transfer rate. This interaction leads to high 
and low states. The mass loss rate of the WD wind (with a 
high velocity of > 1000 km s"') reaches as high as Mwind ~ 
1 X 1O~^M0 yr~', being consistent with the radio observation. 
Thus, the mass transfer rate itself may not be constant but 
vary in time, thus being regarded as a time-averaged rate in 
the present paper. 

From the light-curve fitting Hachisu & Kato also estimated 
the WD mass as Mwd ^ 1-25 Mq and the secondary mass 
to be Mms ~ 3.5 Mq, and concluded that V Sge will ex- 
plode as an SN la in a time scale of ~ 1 x 10^ yr. Since the 
present orbital period of V Sge is 0.51 days dHerbig et al. I 
Il965t [Patterson et al. 1(199 8'). its position in the orbital period 
vs. donor mass plane in Figure [3] indicates C[ > 0. Thus, we 
may regard binaries in the wind phase as "V Sge type stars." 

5.2. Paucity of Progenitor Systems 

The life time of V Sge type stars is typically a few to several 
times 10^ yr, mainly because the time-averaged mass strip- 
ping rate is as high as Mms ~ 1Q~^Mq yr"'. If this channel of 
the WD + MS system produces about four Type la supernovae 
per millennium in our Galaxy (e.g., Cappellaro et al. 1999), 
we should have a chance to ob serve at least sev eral hundred 
V Sge type stars in our Galaxy. ISteiner & Diaz I ([T998i) listed 
four V Sge type stars in our Galaxy and discussed their similar 



properties. Although the masses of the companion stars to the 
WDs are not yet clearly identified, their orbital periods fall in 
the range of 0.2-0.5 days, which is very consistent with the 
orbital periods predicted by our new scenario (see the "final" 
regions in Figs. |5]-[T0]|. However, the total number of V Sge 
type stars is too small (by about two orders of magnitude) to 
be compatible with the new scenario, unless 99% of V Sge 
type stars are hidden. 

The same kind of paucity of the progenitors has been al- 
ready pointed out for supersoft X-ray sources (SSXSs) in our 
Galaxy and is attributed to the Galactic interstellar abs orp- 
tion of supersoft X-rays (iDi Stefano & Rappaportl[l994l) . Di 
Stefano & Rappaport also suggested that circumstellar matter 
m ay play some role in th e obscuration of X-rays. 

iDiaz & Steiner I d 19951) pointed out that soft X-ray flux of 
V Sge is too weak to be compatible with the typical supersoft 
X-ray sources, i.e., at least 2 or 3 orders of magnitude lower 
than that of CAL 87, a prototypical SSXS in the LMC. This 
obscuration may be explained with the absorption of X-ray 
by the stripped matter (or the WD wind itself) and may also 
be related to the observational paucity of the supersoft X-ray 
sources. 

As mentioned in § 1 , iDi Stefano & Kongi (l2003l) reported 
the number of SSXSs in four external galaxies from Chandra 
data. They have estimated at least several hundred SSXSs 
in each galaxy, many of which are obscured by interstellar 
absorption. 

5.3. Angular Momentum Loss by Stripped Matter 

The stripped matter is lost from the binary system with 
some angular momentum. In our treatment, we assume that 
the specific angular momentum (angular momentum per unit 
mass) of the stripped matter is given by equation (|6]l, that is, 
the ablated gas from the companion has the specific angu- 
lar momentum there just at the companion's surface. This 
assumption may be too simplified because the stripped mat- 
ter may get some angular momentum from the binary motion 
during its journey. Here we examine other two cases: one is 
the same as the high velocity WD wind, i.e., 

the other is the slow velocity case, i.e., 

4=1, (20) 

where the stripped ma tter gets large angular momentum from 
the binary torque (see lJahanara et aini2005l for recent three- 
dimensional hydrodynamic calculation). 

For the first case of equation ( fT9] l, we have obtained es- 
sentially the same results as in equation (|6]). If we adopt the 
second case of equation ( |20l ), however, we have common en- 
velope formations in a hundred or thousand years for ci = 3, 
M2,o = 5.0 Mq, and Mwd.o = 1 .0 Mq in Figure|2]regardless of 
Pq. If we start the evolution with ci = 3, M2,o = 4.0 Mq, and 
Mwd.o = 1 .0 Mq, we obtain SN la explosions only for Pq = 2- 
5 days. These results hardly change even if we increase the 
efficiency of mass stripping effect to ci = 10. This is because 
too much angular momentum is removed from the binary for 
the case of equation (l20l i and it makes the separation shrink 
drastically regardless of the ci value. Evolutions with ci = 3, 
M2.0 = 3.5 Mq, and Mwd.o =1-0 Mq result in the same final 
outcome as in equation 

On the other hand, there exist four V Sge type stars wit h 
short orbital periods of 0.2-0.5 days (ISteiner & Diazlll998h . 
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Therefore, we conclude that the angular momentum loss is 
much closer to equation ^ or ( fT9] l rather than equation ( |20] | 
because these V Sge type stars cannot be realized with the 
large angular momentum loss like equation ( l20l i that results 
in formation of a common envelope. 

5.4. Mass Transfer Rate of Simplified Treatment 

Our treatment of thermal time scale mass transfer may be 
too simplified compjired with detailed mass transfer model 
studied by iLanger et al. I (l2000l) and iHan & Podsiadlowski i 
(l2004j) . Han & Pod siadlowski compared o ur results based on 
a simplified model dHachisu et al. Ill999bl) with their detailed 
model calculations, and pointed out that the difference is large 
for lower mass WDs. Although we need detailed mass trans- 
fer model to obtain precise SN la regions, our treatment has 
an advantage of easy and simp le estimation for the SN la pa - 
rameter region. As pointed by IHan & Podsiadlowski I (l2004l) . 
our SN la region thus calculated may deviate from the real- 
istic one for less massive WDs. However, our SN la region 
is probably not so largely different from the realistic one for 
more massive WDs (compare with Fig. 12 of Hachisu et al. 
1999b and Figs. 3 and 5 of Han & Podsiadlowski 2004). 

6. CONCLUDING REMARKS 

Both Cases WIND and CALM originate from the systems 
with massive donors, i.e., young population. It would be 
important to make some comparisons with the observational 
data, such as frequency and population. The red hatched re- 
gions in Figures|5]|7] and|9]indicate a region in which the pro- 
genitor explodes at fjeiay < 100 Myr Also the dashed line and 
the dotted lines correspond to f delay = 200 Myr and 400 Myr, 
respectively. We see in Figure (TTj that Case WIND and thus 
SNe la/Iln (Ila) are realized by the very young system with 
fdeiay< 100-200 Myr 

If MwD.o 0.9 Mq, we have almost no region of Case 
WIND, different from the cases of Mwd.o > 10 Mq. If 
all the WD + MS system with Mi^o > 3 - 6 Mq (ci = 3), 
-^WD,o ^ 1.0 Mq (Mi.o > 6.5 Mq), and Pq ^0.5-2 days pro- 
duces SNe la/IIn (Ila) events (Table|2]i, the frequency of these 
events is estimated to be ~ 5% (including both the WD + MS 
and WD + RG systems with their total number ratio of 4:2). 

A group of Type Iln SNe such as SNe 1997cy and 1999E 
show a ver y similar spectroscopic and ph otometric features to 
SN 2002ic (IWang et al. l2004l:lDeng et a l. 2004; Prieto et al. I 
|2007|) . If these are in fact all Type la/IIn (Ila) SNe, their fre- 
quency can be estimated to be 51} % (Prieto etal. 2007), 
which is consistent with the above estimate. 

Type la supernovae play a key role in astrophysics, and thus 
our progenitor model has important implications. Our model 
depends essentially on the parameter of stripping effect, ci, 
which depends on the properties of WD winds, such as as- 
phericity, velocities, and the efficiency of energy conversion. 
Also we calculate the mass transfer rate using the simple ap- 
proximate binary models. In order to improve these param- 
eterization and approximations, we need multi-dimensional 
hydrodynamical simulations, which are beyond the scope of 
the present study. In the present approach, we constrain the 
ci parameter observationally, and estimate C[ 7- 8 and ci ~ 
1.5-10 from the analysis of V Sge and RX J0513.9-6951, 
respectively. With keeping in mind the necessity of further 
theoretical and observational studies to confirm our new pro- 
genitor systems, we summarize the basic results of our new 
SN la scenario: 

(1) Mass-accreting WDs blow an optically thick wind when 



the mass transfer rate to the WD exceeds the critical rate of 
Mcr 1 X 10"*'Mq yr-^ The WD wind collides with the 
secondary's surface and strips off its surface. If the mass- 
stripping effect is efficient enough, the mass transfer rate to 
the WD is attenuated and the binary can avoid formation of 
a common envelope even for a rather massive secondary. In- 
cluding this mass-stripping effect into our binary evolution 
model of the WD + MS systems, we have found a new evo- 
lutionary scenario, in which a companion as massive as 6- 
7 Mq can produce an SN la for a reasonable strength of mass- 
stripping effect, say C[ ^3. 

(2) We have followed simplified binary evolutions and ob- 
tained the SN la region in the logPo-M2.o (initial orbital pe- 
riod - initial donor mass) plane. The newly obtained SN la 
region extends to massive donor masses up to M2,o ^6-1 Mq 
for Pq ^ 0.5 - 10 days, although its extension depends on the 
strength of mass-stripping effect, ci, i.e., M2.0 ~ 7-8 Mq for 
ci = 10, M2 ~ 5 - 6 Mq for ci = 3, and M2 ~ 4 Mq for 
ci = l. ■ 

(3) We have estimated that the SN la birth rate in our Galaxy 
is t'wD+MS ^ 0.004 yr"' (for ci = 3), which is consistent with 
the observation. The rates of young populations, i.e., fdeiay < 
100 Myr and fdeiay < 200 Myr, ai-e about 50% and 80% of 
the total SN la rate of the WD + MS channel. These short 
delay times of SN la progenitors are consistent with the recent 
observational suggestions that a half of SNe la belong to such 
a very young population as the delay time of f delay ~ 10*^ yr 

(4) Another channel of the WD -1- RG syste m shows a broad 
distrib ution of the delay time over 2-3 Gyr dHachisu et al. I 
Il999ah . thus the two (WD + MS and WD + RG) channels 
yield a bimodality of the delay time distribution. 

(5) The stripped-off material is probably distributed on the 
orbital plane and forms a massive circumbinary torus (or disk) 
around SNe la. Such ckcumstellar matter (CSM) may be con- 
sistent with the observed CSM feature in SN 2006X. When 
SN ejecta strongly interact with massive CSM, it can explain 
the feature of Type la/IIn (Ila) SNe 2002ic and 2006gj. 

(6) Three different environments of SN la explosions can 
be specified by three different states of WDs just at the SN 
la explosion, i.e., the optically thick WD wind phase (Case 
WIND), steady hydrogen burning phase without optically 
thick winds from WDs (Case CALM), and recurrent nova 
phase (Case RN). In Case WIND, SN la ejecta sti-ongly in- 
teract with massive CSM like SNe la/IIn (Ila) 2002ic and 
2005gj because CSM exists near the SN la. The estimated 
rate of Case WIND is ^ 5% of the total SN la rate, being con- 
sistent with the observational estimate. In Cases CALM and 
RN, SNe show a normal SN la feature because the CSM is 
far from the SN but the ejecta may interact with the CSM in 
a much later phase. SN 2006X may be on a border between 
Case WIND and Case CALM. 
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